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First-principles prediction of a decagonal quasicrystal containing Boron 
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We interpret experimentally known B-Mg-Ru crystals as quasicrystal approximants. These ap- 
proximant structures imply a deterministic decoration of tiles by atoms that can be extended 
quasiperiodically. Experimentally observed structural disorder corresponds to phason (tile flip) 
fluctuations. First-principles total energy calculations reveal that many distinct tilings lie close to 
stability at low temperatures. Transfer matrix calculations based on these energies suggest a phase 
transition from a crystalline state at low temperatures to a high temperature state characterized by 
tile fluctuations. We predict B38Mgi7Ru45 forms a decagonal quasicrystal that is metastable at low 
temperatures and may be thermodynamically stable at high temperatures. 
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While evaluating the stability of crystalline phases 
competing with metallic glass formation 0, we discov- 
ered a number of previously unrecognized decagonal qua- 
sicrystal approximants 0- These are ordinary, though 
complex, crystals whose local structural motifs may be 
naturally extended to form a quasiperiodic structure with 
an axis of 10- fold rotational symmetry. We created a se- 
ries of hypothetical quasicrystal approximant structures 
based on these motifs and calculated their total energy. 
Based on our findings we propose that B3§Mgi7Ru45 
should posses a decagonal quasicrystal state that is at 
least metastable, and potentially even thermodynami- 
cally stable at high temperature. 

This prediction is noteworthy because: (1) the com- 
pound contains a substantial amount of Boron, which has 
been speculated to form quasicrystalline structures 0, 
QL IE El though none are yet known; (2) the predicted 
structure is quite different from the established struc- 
tures of Al-rich decagonal quasicrystals yj l§| in that it 
is an intrinsic ternary, while known Al-rich decagonals 
are essentially pseudobinaries; (3) we predict both struc- 
ture and existence from first-principles, while in prior 
work existence was taken as a crucial experimental 
input for structure prediction; (4) first-principles total 
energy calculations confirm the feasibility of entropic sta- 
bilization 0- 

Two experimentally known B-Mg-Ru crystals jToL Hl| , 
B 4 Mg 2 Ru5 (Pearson symbol oP22, space group Pbam) 
and BnMg 5 Rui3 (Pearson symbol oP62, space group 
Pbam), are decagonal approximants. These compounds 
form through solid state transformation at an anneal- 
ing temperature of T a =1323K. Already a relationship to 
quasicrystals is evident in the Fibonacci numbers of Mg 
and Ru atoms in each crystal's stoichiometry. These crys- 
tal structures are illustrated in Fig. Hexagon (H) and 
boat (B) tiles @ are inscribed on these figures to show 
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FIG. 1: top: B 4 Mg 2 Ru 5 (oP22) bottom: BnMg 5 Rui 3 
(oP62). Color coding: Blue=B, Orange=Mg, Green=Ru, 
Empty= Vacancy. Size coding: Large=upper plane, 
small=lower plane. Shading indicates 1/2 atomic occupancy. 
Solid lines outline outline optimal tiling O. Dashed lines in- 
dicate possible phason tiling flips leading to tiling O' . 



how the structures can be interpreted as quasicrystal ap- 
proximants. The tile edge length a q —A.b A is known as 
the quasilattice constant. 

In B4Mg2Rii5 only the H tile appears, and it is deco- 
rated deterministically with every site fully occupied by 
a unique atomic species. The determinism of this dec- 
oration reflects the strong size and chemistry contrast 
between atomic species. That is the sense in which we 
refer to B-Mg-Ru as an intrinsic ternary. Other ternaries 
that form decagonal quasicrystals are pseudobinaries be- 
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FIG. 2: Ideal tiles and decorations. Tiles shown are H, B, 
S and E. Prototiles for the generalized 2-level tiling are 
shown with dashed lines. 



cause at least two of the elements readily substitute for 
each other (e.g. Co and Ni in Al-Co-Ni). 

All atoms lie on two flat layers. The medium sized Ru 
atoms occupy one layer (z=0), while the large and small 
Mg and B atoms occupy the other (z=l/2). We refer to 
this pair of adjacent layers as a slab. Slabs are stacked 
along the c axis with a 3 A periodicity. 

In our deterministic decoration, large Mg atoms oc- 
cupy every tile vertex. Ru atoms occupy two topologi- 
cally distinct sites: tile edge midpoints "Ru Q ", and tile 
interior sites "Ru/j". B atoms form a network of pen- 
tagons and thin rhombi. All B atoms center trigonal 
prisms formed by Ru atoms ^(|- ^ u /3 atoms belong to 
the shared faces of the pairs of trigonal prisms inside the 
tiles, while Ru Q atoms belong to shared edges of trigonal 
prisms in adjacent tiles. Mg atoms serve to cap non- 
shared faces of trigonal prisms. These trigonal prisms 
are common motif in amorphous and related crystalline 
compounds 0,0. 

The experimentally determined atomic positions of 
B4Mg2Ru5 agree very accurately with the deterministic 
decoration of our idealized tiling model. The RMS dis- 
placements of the experimentally determined positions 
from the ideal positions is 0.086 A. The maximum dis- 
placement is 0.14 A found for Ru Q atoms. 

In the experimentally refined BnMg5Rui3 structure 
|10| . neither the atomic occupation nor the decompo- 
sition into H and B tiles is uniquely determined. The 
atomic occupation is not unique because certain B and 
Mg atoms (shaded) have occupancy 0.5. The decomposi- 
tion into H and B tiles is not unique because an HB pair 
can be interchanged by means of a bowtie tile flip Q (see 
dashed lines in Fig. ^| . Tile flips are examples of local- 
ized phason fluctuations ,14] in quasicrystals, but can also 
occur in crystalline approximants to quasicrystals where 
they become discrete configurational degrees of freedom. 
Because the atomic decoration is not deterministic, there 
is no guide as to how to resolve the region covered by an 
HB pair into its separate tiles. 

Occupancy of the partially occupied sites must be 
highly correlated. For example, pairs of 50% occupied 
Mg sites are only 0.5 A apart. A single Mg atom should 
occupy one out of the two sites at any instant in time. In 



addition, we find the nearby partially occupied B sites are 
strongly correlated with the Mg position. When the B 
atom occupies its preferred site relative to the Mg atom, 
the HB pair resolves uniquely into an H tile and a B tile, 
each decorated as shown in Fig. [3] 

Thus we believe the experiment reveals partial occu- 
pancy as a result of tile flip disorder. This disorder could 
be planar in nature, with the HB pairs resolved randomly 
within a single slab that is then stacked with perfect 
3 A periodicity along the c axis. Alternatively (and more 
likely, as we show below) the disorder could be truly three 
dimensional, with an HB pair in a given 3 A slab possi- 
bly resolving oppositely to a BH pair in an adjacent slab, 
introducing a tiling flip between slabs. 

H and B tiles alone can cover the infinite plane 
quasiperiodically. However, other tile types beyond H 
and B are possible (see Fig. To generate the full set 
of tiles we use, start with an H tile and add a bowtie to 
generate a B tile. Adding an additional bowtie to a B 
tile generates either a star (S) tile or else the tile type we 
call E. This process can be continued indefinitely creat- 
ing ever larger tiles, leading to a generalization of the set 
of 2-level tilings All additional tiles in this family 

can be created from the set of prototiles illustrated with 
dashed lines in Fig. [3 

For 10-fold symmetric quasiperiodic tilings, the com- 
position is rigidly fixed by the ideal atomic decorations 
of the tiles. The Boat:Hexagon ratio of r : 1 (here 
t = (v5+ l)/2 is the golden mean) corresponds to frac- 
tions of atomic species as = 1/r 2 = 0.382, XM g = 
l/(r 2 + 2r) = 0.171 and x Ru = t 2 /(t 2 + 2r) = 0.447. 
This ideal composition is about 0.1% richer in B and Ru 
than the composition of the BnMg5Rui3 crystal. Ow- 
ing to the the intrinsic ternary nature of B-Mg-Ru, the 
composition cannot be adjusted to optimize quasicrys- 
tal formation, in contrast to the case of the pseudobina- 
ries, where adjustments of the composition can move the 
Fermi level to a pseudogap, for example. If we wished 
to make comparable adjustments for B-Mg-Ru it would 
be necessary to move to a quaternary (pseudoternary) 
system such as B-(Mg,Zr)-Ru. 

Given our quasicrystal model we study structural sta- 
bility by calculating cohesive energies in the B-Mg-Ru 
ternary system. Our first-principles total energy calcu- 
lations employ electronic density functional theory us- 
ing the plane- wave program vASP 15]. We use the 
GGA approximation with the VASP-supplied PAW po- 
tentials [l|j. Our energy cutoff and k-point densities 
achieve convergence of energy to an accuracy of better 
than 1 meV/atom. For all structures examined we relax 
atomic positions and lattice parameters. We consider all 
known binary and ternary crystal structures |l l| , a large 
number of hypothetical structures drawn from chemically 
similar alloy systems, and 70 different quasicrystal ap- 
proximants. Our methods are described in more detail 
in Ref. [| and quantitative cohesive energy data is avail- 
able on "the WWW [13 ■ 

These calculations exactly reproduce the known binary 
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phase diagrams of B-Mg and Mg-Ru in the sense that 
all known stable phases lie on the convex hull of energy 
versus compositions, and all hypothetical structures lie 
above the convex hull. For B-Ru, all unknown structures 
correctly lie above, and most known structures lie on the 
convex hull. However, two phases, B3RU7 (presumed sta- 
ble) and BgRun (presumed metastable) lie far above the 
convex hull (by 66 and 143 meV/atom, respectively). In- 
troducing vacancies lowers the energy considerably but 
we have not yet found structures whose energy reaches 
the convex hull. These phases require further theoretical 
and experimental study. The ternary diagram is repro- 
duced exactly. 

We resolved the partial occupancy of BnMg5Rui3 in 
many different ways to find the lowest energy. It turns 
out that introduction of tile flips between adjacent slabs 
is favored, lowering the total energy by 8.4 meV/atom. 
The lowest energy structure is obtained by arranging two 
boat and two hexagon tiles in a single slab as shown in 
Fig- H (we call this optimal tiling O), then stacking a 
second slab above that differs by the tile flips outlined 
by dashed lines. The resulting tiling, which we denote 
O' , is equivalent to the starting tiling, but reflected and 
translated. The optimal structure thus exhibits a 6 A pe- 
riodicity perpendicular to the tiling plane, an alternating 
sequence of O and O' with space group Pnma. We assign 
this structure Pearson symbol 0PII6. 

Other tilings exist within the same lattice parameters. 
Indeed, a different arrangement of two boat and two 
hexagon tiles in a single slab, alternating with a partner 
to yield 6 A periodicity and space group Pnma as before, 
has energy just 0.3 meV/atom above the OO' structure. 
We denote this structure as QQ', where Q differs from 
O by just one bowtie flip. 

Starting from ideally decorated tiles we find significant 
atomic relaxation in the tiling flip regions. In particular, 
the Mg atom displaces by over 1 A towards the bowtie 
center (see Fig. ^ bottom). Adjacent Ru Q atoms relax 
towards the now-vacant vertex. The Ru layers become 
slightly non-flat. Our relaxed positions are all within 
0.05 A of the experimentally determined atomic posi- 
tions for BnMg 5 Rui3. 

Our calculated cohesive energies show that the known 
ternary structures are highly stable (enthalpies of for- 
mation are around 350 meV/atom). However, a great 
many hypothetical approximant structures lie quite close 
to the convex hull, starting about 2.5 meV/atom above. 
In other words, we find a cluster of many distinct but 
very nearly degenerate structures in the vicinity of the 
ideal quasicrystal composition. 

To understand why experiments find a disordered 3 A- 
periodic structure for BuMgsRuis, while our calcula- 
tions find an ordered 6 A-periodic minimum energy state, 
we carry out transfer matrix calculations. These model 
the temperature dependence of the thermodynamic en- 
semble of structures that differ by tiling flips. First we 
construct all 2-D tilings a that fit within a given pe- 
riodic boundary condition in the xy-plane. Each tiling 



describes a slab of B-Mg-Ru with height 3 A along the z 
axis, decorated as in Fig. [21 Then we consider combina- 
tions a(3 in which slab (3 is placed above slab a along the 
z axis, for a total height of 6 A. We only consider com- 
binations which differ by disjoint localized tile flips, so 
the tilings match on all vertices neighboring the flipped 
vertex. Specifically, these are the four vertices forming 
the outside corners of the "bowtie" that flips. For each 
combination a(3, we calculate the relaxed total energy 

U a f3- 

We may think of U a p as comprising: the energy of 
slab a; the energy of slab /3; twice the interaction en- 
ergy of slab a with (3. This factor of 2 arises from the 
periodic boundary condition along the z axis. Within 
this energy model, the total energy for a stack of n lay- 
ers (aict2a>3 . . . ,a n ) with periodic boundary conditions 
of height n x 3 A, becomes E aia2 + E a2a3 + • • • + E an0ll , 
where E a p = \U a p. This approximation neglects inter- 
actions of second- and further-neighbor slabs. 

Defining the transfer matrix elements R a p = 
exp (— Ea/s/ksT), the partition function for a stack of 
n layers is Z = tr R n . In the limit of many layers, the 
free energy per layer / = — fc^Tlogp, with p the largest 
eigenvalue of R. Other thermodynamic quantities such 
as internal energy U, entropy S and heat capacity C are 
given by temperature derivatives of /. 

The broad peak in Fig. [3] illustrates the heat capacity 
per atom for the ensemble of structures that match the 
lattice parameters of BnMg5Rui3 in the xy plane. At 
low temperature the structure is locked into the minimum 
energy configuration OO'. It has a vertical periodicity of 
6 A because tiles flip back and forth with perfect regu- 
larity. A small peak around T=500K indicates the onset 
of the alternate QQ' structure. The heat capacity peaks 
around the annealing temperature T a and well below our 
guess at a likely melting temperature T m «1750K. We 
also plot the entropy multiplied by temperature, a mea- 
sure of the free energy reduction due to tiling fluctua- 
tions. At the annealing temperature TS is around 2.6 
meV/atom. At the same temperature, these fluctuations 
increase the internal energy U by around 1.3 meV/atom. 

From the eigenvector of the transfer matrix we deter- 
mine that about 69% of the slabs a occurring in equi- 
librium are of the optimal type a — O or O'. Starting 
from a slab of type O, the next slab above is of type O' 
with probability around 75%. This leads to a persistence 
length for the 6 A periodic OO 1 sequence of 12 slabs, or 
36 A. As a result, there will be no Bragg peak associated 
with 6 A periodicity, but there should be pronounced 
diffuse scattering. 

Does the broad peak in Fig. |3| indicate a genuine phase 
transition? With finite extent in the xy plane, our trans- 
fer matrix describes an effectively one-dimensional sys- 
tem that is incapable of a genuine phase transition. A 
phase transition might exist in the thermodynamic limit 
of infinite extent in the xy plane. This cannot be done 
using ab-initio methods, which are already strained by 
the demands of the 116 atom double slabs of the OO' 



4 



0.02 



0.01 




FIG. 3: Thermodynamic properties of BnMg 5 Rui3 calcu- 
lated from transfer matrix for an infinite stack of single unit 
cells. Left axis, heat capacity C (solid). Right axis, T * S 
(dotted), U (dashed). Arrows mark annealing temperature 
T a and estimated melting temperature T m . 



model. Instead, we approximated the energetics using 
a crude "tile Hamiltonian" [Lsl Il9j that enforces tiling 
constraints and assigns an energy preference for tiling 
flips between adjacent slabs. We evaluate the tile Hamil- 
tonian energies for supercells of the basic BuMgsRuia 
structure and find that the heat capacity peak narrows 
and diverges as system size grows. The existence of a 
genuine phase transition seems clear, although we do not 
know its temperature very accurately owing to our crude 
tile Hamiltonian approximation. 

Could tiling disorder, whose prevalence at high tem- 
perature we have demonstrated, stabilize the quasicrys- 
talline phase entropically 9]? To answer we must com- 
pare the entropy difference between the quasicrystal and 
competing crystal phases with corresponding energy dif- 
ferences. Substituting larger approximants for the true 
quasicrystal, we find the quasicrystal is high in cnerg 
by 2.5 meV/atom. The random tiling hypothesis 
suggests the entropy of the quasicrystal is greater than 



that of small approximants. We already found TS of 
2.6 meV/atom for the BnMg 5 Rui3 phase, but we need 
the entropy difference between that and the quasicrys- 
tal. Since we cannot directly calculate the quasicrystal 
entropy, we take 2.6 meV/atom as indicating the magni- 
tude of the probable difference in TS at temperature T a . 
This rises to 4.0 meV/atom at T m . 

Since both energy and entropy are of comparable mag- 
nitude, we have verified from first-principles the feasibil- 
ity of entropic stabilization. We predict that quasicrys- 
talline B-Mg-Ru is on the verge of stability, so effects left 
out of our considerations could tip the balance. A better 
estimate of the quasicrystal entropy could be achieved 
by finding an accurate tile Hamiltonian then evaluating 
the entropy of this model for large approximants. We 
should include the increase in internal energy U caused 
by tile fluctuations. Further, we should include atomic 
vibrations (phonons) which contribute their own entropy 
and also can be expected to modify the values of terms 
in the tile Hamiltonian. 

In the event that further theoretical (or experimen- 
tal) study shows the B-Mg-Ru decagonal phase is not 
stable at high temperature, moving to a quaternary sys- 
tem (e.g. substituting a small fraction of Mg with Sc 
or Zr) could stabilize it by raising the energy of com- 
peting crystal phases relative to the quasicrystal. Re- 
gardless of the stability at high temperature, we expect 
the ternary quasicrystal to occur metastably at low tem- 
perature, because the energy differences between crys- 
talline and quasicrystalline structures are small, and the 
high entropy of the quasicrystal implies there are a large 
number of structures into which the system could freeze 
out of equilibrium. Thus, we predict the occurrence of 
decagonal quasicrystals in B3gMgi7Ru45. 
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